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Phase-Imprinting of Bose-Einstein Condensates with Rydberg Impurities 
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We show how the phase profile of Bose-Einstein condensates can be engineered through its in¬ 
teraction with localized Rydberg excitations. The interaction is made controllable and long-range 
by off-resonantly coupling the condensate to another Rydberg state with laser light. Our technique 
allows the mapping of entanglement generated in systems of few strongly interacting Rydberg atoms 
onto much larger atom clouds in hybrid setups. As an example we discuss the creation of a spa¬ 
tial mesoscopic superposition state from a bright soliton. Additionally, the phase imprinted onto 
the condensate using the Rydberg excitations is a diagnostic tool for the latter. For example a 
condensate time-of-flight image would permit reconstructing the pattern of an embedded Rydberg 
crystal. 

PACS numbers: 03.75.-b, 03.75.Lm, 32.80.Ee, 32.80.Qk, 


Introduction: The imprinting of tailored phase profiles 
onto the complex order parameter of a Bose-Einstein con¬ 
densate (BEC) [T] is a versatile tool for the creation of 
topological states such as solitons m [3] and vortices [T] 
or even Skyrmions [4]. More generally imprinting allows 
the transfer of the BEC into a desired state of atom flow 
or motion. Typically the phase is generated purely op¬ 
tical using lasers. We propose to engineer phases ex¬ 
ploiting the interactions between BEC atoms that are 
Rydberg dressed [SHIS] and resonantly excited Rydberg 
atoms pAlflS] . This phase imprinting will rely on the 
matter-wave coherence of the condensate, and thus rep¬ 
resent an instance of genuine Rydberg-BEC physics. 

We show that phase imprinting creates a versatile in¬ 
terface between ultra-cold Rydberg- and BEC physics. 
Firstly, it allows entangled Rydberg states [T9]-[2T] to be 
mapped onto the many-body wave function of the con¬ 
densate. As one example we discuss how to turn an 
atomic Bell state | + ) = (| Rg ) + | gR )) / (with two 
atomic electronic states \g), |R)) into a spatial meso¬ 
scopic superposition state in the position of a single BEC 
bright soliton, akin to the proposal of [22]. Secondly, 
phase imprinting represents a tool to probe Rydberg elec¬ 
tronic states via their effect on a condensate diiMiT]. 
To demonstrate, we show signatures of Rydberg-crystals 
[28] in expected condensate time of flight spectra. 

Our scheme relies on weakly admixing Rydberg char¬ 
acter to all the atoms in a condensate cloud through far 
off-resonant laser coupling between their stable ground 
state and a highly excited Rydberg state 00 , as ex¬ 
perimentally demonstrated for two atoms m- All these 
’’dressed” atoms then interact with some previously pre¬ 
pared, fully excited atoms in a different Rydberg state, 
referred to as impurity or control atom depending on 
whether they are inside or outside the condensate. The 
atomic species of impurities may be identical to conden¬ 
sate atoms. 

The interaction involving impurities can be made dom¬ 


inant over the simultaneously induced long-range conden¬ 
sate self-interaction 00. Using this we demonstrate the 
imprinting of sizeable phases for realistic parameters. 
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FIG. 1: (color online) (a) The setup. Several Rydberg ex¬ 
cited impurities (red balls) are embedded in a Bose-Einstein 
condensate. The condensed atoms (green balls) mainly in¬ 
teract with impurities through dressing induced long-range 
interactions of characteristic range rc (red dashed circles). At 
distances d within the Rydberg orbit, d < ro, collisions with 
Rydberg electrons are also relevant (blue dashed circles). Af¬ 
ter an imprinting period, the condensate phase (blue shades) 
will be modified only in the vicinity of impurities (dark blue). 
(b) Level scheme for three representative atoms. Interactions 
Vrr give rise to long-range interactions between dressed con¬ 
densate atoms as in [iH. Interactions VrR and VgR give rise 
to stronger potentials between dressed condensate atoms and 
impurities, which are the main focus here, (c) Sketch of these 
potentials U^^ (red, see text) and VgR (blue) near one of the 
impurities, a = U/(2A) quantifies the degree of Rydberg 
admixture. 


The condensate phase is also affected by direct colli¬ 
sions of the Rydberg electron with condensate atoms in 
the ground state, without dressing coupling. This can be 
exploited to visualise the electron orbital through con- 
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densate densities [H HI]. In contrast the imprinting 
through dressed interactions discussed here extends the 
spatial scale, smoothness and controllability of phase pro¬ 
files. The use of Rydberg impurities and effectively one- 
or two dimensional (ID,2D) condensates in this article 
circumvents some of the interaction strength and many- 
body related problems of dressing discussed in m- 
Interactions between Rydberg impurities and dressed 
atoms: Consider a gas of N Rb atoms with mass M 
at locations R^, as depicted in Fig. (a). They may be 
in either of three electronic states: ground state | ^), Ry¬ 
dberg state |r) = \ns) or Rydberg state |R) = \n's), 
with principal quantum numbers z/, z/', where v < v' 
as shown in Fig. (b). The label s implies angu¬ 
lar momentum I = 0. We write the system Hamil¬ 
tonian Htot = Hq Hdress + ^int: usiug the notation 

= \k){k' \ ^ where /c, k' G {^,r, R} and acts on 
atom n only: 

N 

Ho = Y.[- + W^(Rn)<7W] , 

n=l 

N 

Rdress = ^ + h.C. - , 

n=l 

N 

Hint= E (1) 

a,bE{g,rR} n^m=l 

Ground state atoms experience an external trapping po¬ 
tential IF(Rn). The state | 5 ^) is coherently coupled to 
|r) with Rabi frequency Il{t) and detuning A. Defin¬ 
ing dnm = iRn ~ Rm 1 5 wc take vau-dcr-Waals inter¬ 
actions between two Rydberg atoms in states a, 6 as 
Vab{d) = /[{I + Sab)d^] for simplicity, where 

is the dispersion coefficient and Sab Kronecker’s delta. 
Between ground state atoms we assume the usual con¬ 
tact interaction Vgg{d) = gS{d), where g = Anfisas/M 
with atom-atom s-wave scattering length a^. Finally, 
ground state and Rydberg electrons interact via Fermi 
pseudo potentials Vg^r/nid) = with Vq = 

2'Kl?aajma [29], electron mass rUe, electron-atom scat¬ 
tering length Ue m and Rydberg orbital wave func¬ 
tion Our examples will be based on Rb with 

Oe = —0.849 nm, assuming Rydberg states |r) = |555') 
and I R) = I 76S ) [3T] . 

Now consider a scenario where Aimp of the N atoms 
have been excited to the impurity Rydberg state | R), de¬ 
noting their locations by {x4 C {R4. Many different 
random or deterministic location patterns can be created, 
depending on the method of excitation. We do not con¬ 
sider the excitation step, but refer to the literature on se¬ 
lective optical access [321134] , the exploitation of blockade 
effects [28j 1351143] or on condensate density dependent 
energy shifts [T6l[23]. Interactions Vrr between two Ry¬ 
dberg impurities are important in the stage of impurity 


placement but can subsequently be neglected in the ex¬ 
amples discussed here. The light-atom coupling in ^dress 
will cause long-range interactions for all atoms, which 
would otherwise be present only among Rydberg excited 
atoms (I r ), | R)). Assuming far off-resonant coupling 
between | g) and | r ), so that \a\ ^ 1 for a = Cl{t )/(2A), 
we determine these interactions using fourth order per¬ 
turbation theory in ft{t). Calculating the energy shift 
AE'gR of the state | 0) = | gR) where the atoms 

at locations are in | R) and the other N = N — Aimp 
atoms in I ^ ), we obtain AE^n = a^E(^ 2 ) +<^"^^( 4 ): where 

^( 2 ) = AEf(l-Em"P Vi?(|Rn-x„|)/A)-i. We find 
^( 4 ) in comparison negligible, see [44]. 

In this article, impurities will only affect dynamics for 
very short times, such that their motion can be assumed 
frozen in space [45], and they also do not undergo state 
changes. Similar to [6] we merge the effective interac¬ 
tions obtained through the laser dressing with the direct 
interactions between atoms contained in 0 i4g) 

to arrive at the following effective Gross-Pitaevskii equa¬ 
tion (GPE) for the dressed and condensed ground-state 
atoms in the presence of Rydberg impurities: 


= ( - 4-v^ + w(R) +5|4R)I^ 


dt 


2m 


+ 


N-„ 


f/i2)(R,{x„},f)+ ^ yo|^-(|R-x„|)|2 j4R), 

N- 

- Y, K-k(|R - x„|)/A)-i. (2) 


Here the presence of a few impurity atoms that are 
fully in a Rydberg state causes strong, long-range inter¬ 
actions with the remaining atoms, which can be treated 
as external single body potential for the condensate. The 
corresponding terms in 0 are in square brackets and will 
be the central tool of the present work. Note that the 
dominant part E(^ 2 ) is orders of magnitude larger than the 
dressed interaction between condensate atoms, causing 
quite different physics than the latter [6]. The induced 
potentials are sketched in Fig. [^(c) as red line , us¬ 
ing signs VrR / A < 0) and blue line (| T (| R|) p) for a single 
impurity. Either potential is associated with an impor¬ 
tant length scale. The plateau of the dressing induced po¬ 
tential extends to the critical radius Tc = ^ 

which also sets the width of the region of potential drop. 
The extent of the direct interaction potential VgR is the 
radius of the Rydberg electron orbital of the impu¬ 
rity, ro ~ with ap the Bohr radius. We focus on 

parameters for which molecular resonances are avoided 
and also Tc > rp [46H48]. Although included in our so¬ 
lutions of Eq. Q [49], the direct interactions VgR then 
play a minor role. It has been shown that many-body 
perturbative calculations as used here are valid only as 
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long as there would be much less than one Rydberg ex¬ 
citation blockade sphere (A^bi ^ 1) [HllSn], which will be 
satisfied here. 

Our applications of Eq. il) to phase-imprinting in¬ 
volve two dynamical stages: In a first short stage of 
duration rimp ^ 10/is, the condensate order parameter 
(j) acquires a dynamical phase such that 0(R) ^ 

exp [i(^(R)](/)(R) with (^(R) = -J7g^g^(R, {xm})Timp. The 
time Timp and strength of are such that other en¬ 
ergies can be neglected. Only in this stage are dressed 
interactions enabled though ^dress- In a much longer sec¬ 
ond stage {t ^ 10ms), the condensate evolves according 
to the usual GPE, and the initially imprinted phase pro¬ 
file is typically transformed into condensate flow and/or 
density variations. We consider two examples that high¬ 
light the main strengths of Rydberg phase imprinting: 
Transferring entanglement from a Rydberg system onto 
a BEG, and inferring the geometry of a collection of Ry¬ 
dberg impurities in a cold gas. 

Entanglement transfer: We first consider a ID arrange¬ 
ment of a ^^Rb BEG bright soliton (see Ref. [5TJ [52] 
and references therein) with N = 400 atoms, located 
between two individual atoms, which are each tightly 
trapped in their own optical tweezer at xi ^2 with po¬ 
sition spread aeon = 0-05 jam. The atoms outside the 
condensate are referred to as control atoms. As shown 
in Eig. the control atoms are separated by a distance 
D = \xi — X 2 \ = 3 fim. The soliton requires attractive 
contact interactions, ^ < 0, which can be achieved using 
the Eeshbach resonance [53] at 5 ^ 155 G in ^^Rb [54] . 
we assume as{B) = —5.33 nm [55]. Only the control 
atoms are now driven into the Rydberg state | R) = \u') 
under blockade conditions, resulting in an entangled two- 
body state 1+) = {\gR) + \ Rg))/V2. Subsequently 
we enable the dressing coupling Q{t) to the state |r) 
for the bulk soliton, resulting in the dressed potential 
sketched blue in Eig. (a,b), which depends on whether 
the left or right control atom was originally excited. After 
an adiabatically enabled and disabled imprinting period 
^imp = 36/is, with Q{t)/h = 3 MHz, A/h = —500 MHz, 
the condensate has acquired the phase profile shown in 
red. Hollowing de-excitation of the control atoms, we al¬ 
low Tmov = 2 ms of free evolution according to the first 
line of Eq. (HI, i.e. Q{t) = 0. After Tmov, the soliton has 
moved by about 2 fim to the left or right depending on the 
imprinted phase profile, as shown in Eig. [^(c,d). Let us 
denote the many-body wave function of the gas for these 
two cases as | T )ieft/righf entire process should be 
quantum coherent, since the initial imprinting happens 
well before a control Rydberg state or dressed Rydberg 
state from the condensate would decay [56], resulting in 
a final many-body state l«'+) = (l«'>left + l^)right)/v^- 

The tightly trapped control atoms m and the small 
length scales [32] are technical challenges for the above 
proposal. However, when these are overcome, one obtains 


a mesoscopic entangled state [22] [581163] , where the entire 
soliton of N atoms is in a superposition of two different 
locations. As proposed in [22l[64] . 

The superposition nature of the resulting state can be 
proven interferometrically [22] [64] [65] upon recombina¬ 
tion, for which one would additionally place the soliton 
into a weak harmonic trap 1T(R). De-coherence pro¬ 
cesses during the creation of such a highly entangled 
many-body state limit TV, but are small for our choice 
here [66]. A full quantum many-body treatment includ¬ 
ing coherence between control atom states and conden¬ 
sate atoms may be subject of further research. 



FIG. 2: (color online) Entanglement transfer from two control 
Rydberg atoms onto a mesoscopic BEG cloud. The position 
space density of control atoms is shown as shaded curves, blue 
for state | R) and gray for state \ g). The condensate density 
for a soliton is shown as thick black line, (blue) Dressing 
potential (red) condensate phase (f after imprinting. 

We show all quantities with arbitrary normalization to fit the 
same axis. (a,b) Initial state at t = Timp, (c,d) final state at 
t = Timp + Tmov- (a,c) show the control atom configuration 
\Rg), (b,d) show \gR). For each configuration, we model 
condensate evolution separately using Eq. ([^ [55]. The dot¬ 
ted line in (c,d) is for shifted control atom positions Xi-\-2acon- 

Rydberg erystal imprinting: The maximally entangled 
state is but one example of entanglement arising 

due to strong Rydberg-Rydberg interactions. Another 
example is given by spatially ordered (crystal) structures 
formed by large numbers of Rydberg excitations in a cold 
gas [28] [40] [42] [67I[69] . We show now that Rydberg- 
phase-imprinting in the presence of such structures leads 
to condensate momentum spectra that allow the recon¬ 
struction of the locations of impurities. 

Let us consider a 2D model [70] of a BEG confined 
in a pancake shaped harmonic trap IT(R) = m[uj‘^{x‘^ + 
^^)+Co’^ 2:^]/2, with frequencies uJz ^ and = 5.5 nm, 
thus ^ > 0. Using a scheme as discussed in [28], A^imp 
impurities can be arranged, for example, in a crystal like 
structure within the condensate cloud. 

Eor a distribution of impurities as shown in Eig.[^(a), 
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FIG. 3: (color online) Determining the spatial arrangement of 
Rydberg impurities (shown as blue crosses) via phase imprint¬ 
ing in a 2D ®^Rb BEG of N = 100 atoms in a trap with Ur = 
(27r) 2Hz and Uz = (27r) lOOHz. Golorbar for (a,c,d) see (d), 
where no is the respective peak density, colorbar for (b,e), see 
(e). (a) Initial condensate density p — |0(R)|^. (b) Gonden- 
sate phase p following phase imprinting {t — Timp). (c) Gon- 
densate density shortly after imprinting (t = rimp + 12ms). 
(d) Final momentum spectrum |0(/c)| (t = 60 ms), (e) Phase 
prohle reconstructed from (d) as described in the text. 


we numerically solve Eq. (I^ [TTJ [72] enabling the po- 
tentials and VgR |l9] for a short imprinting period 
Rmp = 18.5/is only, using Q/h = 4 MHz, A/h = —150 
MHz. This is followed by evolution under the influence of 
the contact interactions, but with disabled harmonic trap 
and the impurities assumed removed via field-ionisation 
[ElEs]. At some final time where momentum spectra no 
longer significantly change, we plot the expected time-of- 
fiight (TOE) images in Fig. We also show the position 
space density shortly after phase-imprinting. 

If the effect of atomic collisions, described by the non¬ 
linear term ^|0(R)p, is not too large, the final momen¬ 
tum spectrum is roughly the Fourier transform of 0(R) = 
Y^nini(R) exp[i(^(R)], where nini(R) is the known initial 
atom density in the trap, and (p{Il) the phase profile gen¬ 
erated through imprinting and shown in Fig. |(b). A 
standard phase-retrieval algorithm HIITSI is then able 
to recover the phase profile as shown in Fig. (e) from 
which impurity positions can clearly be inferred. The 
algorithm relies on iterative Fourier-transforms involv¬ 
ing two known quantities: the final time-of-fiight image 
from which the modulus of the condensate order parame¬ 
ter is extracted and the initial condensate density 

^ini(R)- We briefly describe it in [44]. 

We find that simple phase retrieval fails for the case of 
Fig-i but larger condensate densities due to condensate 
self-interactions. This might be remedied by more so¬ 
phisticated variants of the phase retrieval algorithm m, 
or modifications of atomic interactions using a Feshbach 


resonance [53]. One could then use larger atom clouds, 
in which case impurity locations can be recovered from a 
single image as in Fig. (d) without the need for image 
alignment in an ensemble average, as in [42]. 

The setup just discussed complements Rydberg crystal 
detection based on single-atom addressing [42] or electro- 
magnetically induced transparency [24f[26] by working 
with a bulk gas and moving the signal from the light to 
the atomic density. Beyond crystal detection, it enables 
phase-profiles that are otherwise difficult to achieve, for 
example those akin to Fig.[^(b) arising from a crystalline 
impurity distribution on the surface of a 3D sphere. 

Conclusions and outlook: We proposed a novel phase 
imprinting technique for Bose-Einstein condensates, 
employing long-range interactions between condensate 
atoms and embedded Rydberg excited impurity atoms, 
created by coupling condensate atoms far off-resonantly 
to another Rydberg state. The scheme offers function¬ 
alities beyond existing imprinting methods, as it allows 
mapping of entanglement from few-body Rydberg states 
onto the whole atom cloud, and strengthens BEC as a 
diagnostic tool for detecting Rydberg excitations in an 
ensemble of atoms. We illustrate the former through a 
proposal for the creation of a mesoscopic entangled state 
in the position of a cloud of atoms and the latter by ex¬ 
ploring the link between Rydberg crystal structures in a 
condensate, and momentum space spectra after phase¬ 
imprinting. 

Combing the techniques discussed here with imprint¬ 
ing effects by a Rydberg electron in a larger orbital may 
offer additional possibilites, due to the unusual shape of 
the Rydberg orbital pTUIS] . Other interesting physics 
might arise from the interplay of phase-imprinting and 
controlled impurity motion [76U8Q] . 

We gladly acknowledge fruitful discussions with Igor 
Lesanovsky, Thomas Pohl, Rick van Bijnen and Shannon 
Whitlock, as well as EU financial support received from 
the Marie Curie Initial Training Network (ITN) COHER¬ 
ENCE”. W. L. is supported through the Nottingham 
Research Fellowship by the University of Nottingham. 


Supplement: Phase-Imprinting of Bose-Einstein 
Condensates with Rydberg Impurities 

Dressed ground state - impurity interaction: Here we 
analyze effective interactions arising from VrR and Vrr 
in (1) of the main article perturbatively. Consider 
H' = h.c.]/2 the perturbation and 

Hq = ^tot ~ H' the unperturbed Hamiltonian. We define 
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states 


|0) = \9l 92-■■ 9 nRi ■ ■ ■ -^Afimp) ’ 

N) = \9i > 

\ij) = \9i 92---ri...rj ... • • • -RjVimp) ’ (3) 

with zero-, one- and two atoms in the Rydberg state | r), 
and a fixed number of A^imp impurities in the Rydberg 
state I R). 

Let us introduce the shorthand = T46(|I^n —R/l)* 
We then calculate corrections Eq to the energy of the 


state I 0). Up to fourth order perturbation theory in H' 
we obtain Eq = + o?E(^ 2 ) + <^^^( 4)5 where we define 

the unperturbed energy = 0. We find to leading 
order 


N 


a‘ 


'■Ei2)=a^^Y. 


1 - VA 


(4) 


Since the dressing parameter a = U/( 2 A) ^ 1 , this 
will be the dominant consequence of laser dressing. The 
fourth order correction is 




N 

E 

n,n';n7^n' 


(A - - E^rf f^;/L)(2A - - Kn') 


(A - E^=f + KD - Kn') 


N 


-E 


E (A - urR){A - j:lTr Ktr 


(5) 


To understand the physics contained in £^( 4 ) and ulti¬ 
mately find a simple approximation, we analyse £^( 4 ) in 
some limiting cases. When all interactions vanish: 

lim = -Aa^N, ( 6 ) 

as expected for the fourth order light shift of N non¬ 
interacting background atoms. 

If we let all background atoms interact, but remove the 
effect of impurity atoms (e.g. placing them far away), we 
obtain 


This is the fourth order contribution to the effective inter¬ 
action of one background atom and the impurity atom, 
the leading contribution of which is a‘^E(^ 2 )- 

Beyond these limits, ^ must contain nontrivial three- 
body physics, as becomes clear when we consider the 
impact of an impurity atom on the effective interaction 
between two background atoms, 0 : Whenever either of 
two atoms 1 or 2 separately falls into the blockade radius 
of the impurity, it cannot reach its Rydberg state | r), 
and the interaction 0 must be suppressed. 


N 


lim a^Eu^ = -Aa^N + Aa^ 


E 

n^n' 


Tjrr 
^ nn' 


2 A - 

^ nn' 


(7) 


This is simply the interaction between two Rydberg- 
dressed ground state atoms that is known from [ 6 ]. 

If instead we neglect interactions between background 
atoms the expression becomes 


lim a^L^( 4 ) = 


r{rr) _ 


-Aa^^ 




(1 - U^^r^^/Af 


( 8 ) 


Here we only wish to consider the dominant conse¬ 
quences of the induced interactions. The interaction 0 
is suppressed by compared to Q , however due to the 
much larger number of background atoms compared to 
impurity atoms, care has to be taken when estimating 
the importance of the term. 


In order to simplify this estimate and (for the sake of 
completeness) allow an inclusion of the effects of 0 in 
a Gross-Pitaevskii description, we propose the following 
simplification of £^( 4 ), based on the discussion above 


N 

«"^( 4 )=E 




N 


(1 - Em™P Va)^ 


Aa^ E 


Tjrr 
^ nn' 


2 A - 

n,n';n'7^n 


(1 - E™™’’ u ^^^/ a )' 


n 


(9) 
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The term in square brackets is a phenomenological 
screening factor that suppresses effective background- 
background interactions whenever atom n is too close to 
an impurity (we see shortly why we do not include atom 
n' in the screening factor). It can be easily verified that 
(§ agrees with (§ in the limiting cases discussed earlier. 
Since these span most of parameter space, the approxi¬ 
mation of the original ^ works usually well. We con¬ 
sider deviations for a special case with two background 
atoms and one impurity atom. Deviations are expected 
for geometries where the distances of all three atoms are 
comparable and of the order of the blockade radius. This 
is confirmed in Fig. Deviations are small enough for 
our purpose, which is merely to ascertain that fourth or¬ 
der interactions © can be neglected for parameters here. 



FIG. 4: (color online) Benchmark of our simplification for the 
fourth order dressed potential (a) Full potential 

(a^F( 4 ) ^ for fixed d 2 R = 4.7/im. (b) Relative difference 

SE = |F(4) — F(4)|/F(4). 


In fact we only employed Eq. © to confirm that for 
all cases discussed in detail in this article, can 

be neglected relative to (a^E( 2 ). The continuous space 
version of (§ would give rise to a term 

+ J ( 10 ) 


on the rhs of the GPE (2) of the main article. Here it 
becomes clear why the approximation involved in Eq. 
is advantageous: We can now write the non-local inter¬ 
actions in Eq. (10) as 


/ 




/ 




d^K'[U^J>(R - R0l<^(R0P /scr(R, 

( 12 ) 


where the term in square brackets can be efficiently evalu¬ 
ated via a convolution as usual. Had the screening factor 
in contained n' this would not have been possible. 

Eor all results in the main article, an interaction 
Eq. ( [T^ would give rise to a dynamical phase of ip < 
O.OOStt during rimp- 

Phase recovery algorithm by Gerchberg and Saxton: We 
briefly summarize here the algorithm due to Gerchberg 
and Saxton m for the recovery of phase information of a 
function from just the intensity (density) of the Eourier 
transform of the function. The algorithm and variants 
ca are widely applied, e.g. in x-ray crystallography. 

Consider a pair of function 0(R) and its Eourier trans¬ 
form 0(k) defined by: 

(^(k) = = (27r)-3/2 J (13) 


Eor our application here, 0(R) and 0(k) are the conden¬ 
sate wave function in the position- and momentum rep¬ 
resentation (respectively), but the algorithm described is 
more general. Let us define the modulus and phase of a 
function f = p[f] exp {iip[f]]) with p, (/? G M. 

The problem is to recover (p[0] when only p[0] and p[0] 
are known. This is solved by repeated application of the 
following steps 


fk=^-^[9k], (14) 

f'k = (15) 

5fc=-^[/fe], (16) 

flfc+i = p[^] exp (Mg'k]), (17) 


with 


AT- 


f/i^“)(R, {x^}) = -Aa\l - ^ [/(>■«) (R - x™)/A) 


imp 


-3 


u: 


i^'’)(R, R', {x„}) = C/^7\R - R')/scr(R, {x™}), 


f/i7)(R-R') = Aa' 


tPG/ 

iff ' 

U^^(R - R') 

2A - [/’’'■(R - R') ’ 


iV;. 


/scr(R,{x„}) = (l- ^ (R-X„)/A)-3, (11) 


where k is the iteration number and the algorithm is ini¬ 
tialized with go = p[0]. In Eq. (15) and Eq. (17) we are 
enforcing the known position and Eourier-space density 
constraints. The algorithm usually converges but possi¬ 
bly very slowly, hence improvements exist m- Eor our 
purposes (Eig. 3 of the main article), between 2 x 10^ and 
2 X 10^ iterations of the basic algorithm were sufficient. 
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